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We present an expression for the transition probability between Dirac or Majorana 
neutrino flavors obtained from first principles within quantum field theory. Our 
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derivation is based on a standard quantum mechanical setup and includes the specific 
mechanism of neutrino production only in as much as it specifies the initial state. Our 
expression for the transition probability reproduces the usual formula in the plane 



also allows a simple understanding of the decoherence of the oscillations and of the 
question of the arrival times of the different neutrino mass eigenstates. We show 



numerical examples for the case of two neutrino generations. 
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INTRODUCTION 



00 

Neutrinos are hoped to enlighten the search for new physics at both the cosmological and 

£> ' I I 

.£h ! microscopic level. The discovery of neutrino oscillations at the SuperKamiokande [lj and 

SNO j^, 3] detectors prompted a detailed examination of the theory of neutrino oscillations 
(some early references are Pontecorvo 4|, Kayser Bilenky jf|; see Giunti and Kim [?| for 
further references). 

Some papers were dedicated to obtaining the formula for the oscillation probability from 
first principles. There are some aspects in these derivations that we do not find very precise 
and that we want to address in this paper. 



Beuthe js] and Giunti [9] stress the fact that the pattern of neutrino oscillations depends 
on the mechanism of production and detection, in other words, the oscillation amplitude 
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is different if neutrinos are produced in pion decay or from a charged current interaction. 
We believe that, following standard rules of quantum mechanics, the pattern of oscillations 
should only depend on the initial state of the system. The particular production mechanism 
needs to be specified only to determine the initial state. Analogously, for detection, quantum 
mechanics can provide the probability of detecting a neutrino in any specific final state. In 
other words, the theoretical description of the phenomenon of oscillations should involve 
the initial state, its evolution and a final state representing the detection process. This is 
because of the quantum mechanical properties of the system, which state that different ways 
of producing the same initial wave packet should bring to the same result. 

In their articles, Beuthe [8] and Giunti [0] average their final formulas for the oscillation 
probability over time, justifying this because the arrival time is not detectable. We will not 
perform this average because, from a fully theoretical treatment, oscillations occur both in 
space and time on the same footage. We will show how this feature of neutrino oscillations 
is related to the covariance of the oscillation formula and to the problem of defining an equal 
energy or equal momentum condition for the various mass eigenstates. 



Snellman and Pallin 



10j derived the expression for the oscillation amplitude by prop- 



agating a second-quantized neutrino wave packet in vacuum and in matter, intentionally 
neglecting any theory concerning production and detection of neutrinos. However in their 
paper there is a misuse of flavor and mass eigenstates. In fact, for the construction of the 
wave packets, the authors in [h]] use the neutrino flavor representation and introduce con- 
cepts such as the energy E a of a neutrino with flavor a. Since neutrino flavor states are 
not eigenstates of the free Hamiltonian, the flavor-eigenstate energy E a is not well defined. 
Moreover, the authors in 



lOj use the formula E a = ^p 2 + m^,, where they state that 



is the mass of the flavor a. This is misleading because a flavor eigenstate, being a linear 
combination of different mass eigenstates, does not have a well-defined mass. 

In this paper, we recover the formula for neutrino oscillations in flat space-time using the 
propagation of second quantized Gaussian wave packets. We choose wave packets instead of 
monochromatic plane waves because we want to represent neutrino states that are "localized" 
near the production and later near the detection point, but plane waves are "everywhere" 
at the same time. In addition, wave packets are normalizable while plane waves are not. 
Another reason is that by using them we can study the coherence between eigenfunctions 
with different mass (and thus different group velocities). We also remedy the misuse of 
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flavor states and flavor properties in [h]] by deriving the neutrino oscillation formula using 
mass eigenstates only. After all, it is the mass eigenstates that appear in the perturbative 
expansions of the scattering processes, for example in the Green functions and the Feynman 
rules. 

Finally, we show numerical results for the case of two neutrino generations, for realistic 
mass differences and for the pedagogical case of a large mass difference. The latter case 
illustrates how decoherence suppresses the oscillation phenomenon. 



II. FORMALISM 



Before recovering the transition probability formula we first establish our notation by 
reviewing the description of the Dirac and Majorana free fields in Quantum Field Theory 
in the Heisenberg picture. We also review the wave packet formalism to be used later on. 

A Dirac field $i(x) describing a particle of mass m,- is fll] 

U*) = E / d A&(P. C)«(p, C)e~ ip - X + d\(p, (Hp, C)e i n , (i) 



where 



is a Lorentz-invariant measure and p ■ x = p^x^ — Et — p • x. 

Eq. is a solution of the free Dirac equation [i^^d^ — mj^i) = 0. The operator 
S|(p, C) creates a neutrino of mass m», momentum p and spin polarization £, while d\(p, Q 
creates an antineutrino of mass m,, momentum p and spin polarization (. The operators 
bi(p, C) and di(p, () annihilate the respective particles. The particle states belong to the 
Hilbert space TC, which can be taken to be the Fock space constructed on the vacuum state 
|0) defined by the condition 6«(p, C)|0) = di(p, ()\0) = 0. The spinors u(p, Q and u(p, () are 
defined by 

u(p,C) = VE l (p)+m, l \ X(P ' ], (3) 




(P,C) = y/EM+mt ( ^+- X(P ' I . (4) 
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Here x(P? C) is a two-component spinor of momentum p and fixed polarization ( which we 
take to be independent of m 8 and p. The spinor x(p> is normalized according to 



From Eqs. ([3] -[5]) we get 



xHp,Ox(p,v) = <W ( 5 ) 

«(p 1 C)7MP.l) = 2p% 1 (6) 
«(P,C)7^(P^)=2^, (7) 

«(p,ChV-p^) = o. (8) 

The fundamental anticommutation relations between the fields ipi(x,t) and their canon- 
ical momenta f[;(x, t) = i^|(x, i) are: 

{^(x,t),fl,(y,t)} = {^(x,t),#j(y,t)} = i5 3 {^-y)5 ij , (9) 

{^(x,t),^-(y,t)} = , {n i (x,*),n i (y,t)} = 0. (10) 

From these equations we derive the anticommutation relations between the operators of 
creation b\(p,() and annihilation 6j(p, C) 

{6i(p, 0, St(q, /?)} = (2tt) 3 2^(p)5 3 ( P - q)5 c „<% , (11) 

{S i (p > C),6i(q,»?)} = 0, {St(p,C),6](q,7?)} = 0. (12) 

Similar relations with 6 replaced by d hold for the antiparticle creation and annihilation 
operators df(p, C) and c?i(p, £)■ 

A Majorana field is described by Eq. (DO) with the additional relations 

ipi(x) = Cipi T (x) , di(p,() = 6i(p,C) > ( 13 ) 

where C = ry°7 2 is the charge conjugation matrix. To facilitate the transition to a Majorana 
field, we choose the phases of u(p, £) and v(p, () so that v(p, C) equals the charge conjugate 
of w(p, (), namely v(p, () = Cu T (p, (). The quantized Majorana field is thus explicitly 

= E / ^[^(P ) C)^(p,C)e-^ + &l( P ,C)Cn T ( P ,C)e^] . (14) 
The normalization of the 6j(p, C) an d S|(p, £) operators is the same as in Eq. ffTTj) . 
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We now discuss the states of the system. In the Heisenberg picture free states do not 



A free state |p, £, i) of momentum p, polarization ( and mass rrii is given 



evolve in time 
by 

\pX,i) = -M=L=\0) . (15) 
The factor in the denominator comes from our choice of normalization condition 

{p,C,i\q.,vJ) = £ 3 (p-q)<*cA/ » ( 16 ) 

which we use for both Dirac and Majorana particles. 
The states |p, satisfy the eigenvalue equations 

P»\p,(,i)=P> M \p,(,i) , (17) 
#|p,C,*> = |P,C,*) , (18) 
^n|p,C,*)=C|p,C,*> , (19) 

where the four-momentum operators P M , the number operator N, and the spin operator 
projected onto the polarization axis S n are respectively 

= E / d pp" K(p» 06i(p, o + 4(p, okp, 0] , (20) 

c J 

n = J2 [ d p Kfo c)6i(p, o + <2(p, o*(p, 0] , (21) 

Sn = J2 [ d P{®(Pi(MP,0 + dl{p,()di{P,0], (22) 

if the states are described by the Dirac theory, or 

ptl = J2 [ #P"&ftp»C)fc(p,C) , (23) 
C J 

N = J2 [ dpb\(p, C)Up,C), (24) 
Sn = Yl [ dp(bt(pXMP,0, (25) 

in the Majorana theory. 

As a consequence of Eq. ( IT6l) . the states |p, £, «) belong to a continuum spectrum and are 
not normalizable. Following standard procedure, we introduce wave packets in Fock space. 
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Let 

l c i> = S / rf3 P a (P;C)|p,C,*>- (26) 

This is a wave packet specified by the function a(p, £), which represents the probability 
amplitude in momentum space for the state Notice that wave packets in the Heisenberg 
representation do not depend on time. The wave packet is normalized according to 

(ci\ci) = J2 [ d 3 p\a(p,C)\ 2 = l- (27) 

Another way to see that a(p, () is the momentum probability amplitude is the expression 
of the momentum expectation value 

(c\P»\c) = J2 j rf 3 p^|a(p,C)| 2 - (28) 
C J 

This expression applies for both Dirac and Majorana particles. 

We will discuss the use of the momentum-space distributions a(p, Q to derive the prob- 
ability amplitude in Section HVl 

The first-quantized wave function corresponding to the free wave packet |cj) in Eq. 
can be computed from Eq. (pQ) and (TT5T) to be 

The wave function ^(x, t) satisfies the Dirac equation for any choice of the momentum-space 
probability amplitude a(p, C). The relativistic wave function ^(x, t) is not the probab ility 
amplitude in position space. The latter, in the interpretation of Newton and Wigner 
is given by the Fourier transform of a(p, £), 

5(*>0 = / 7^«(P>0e ip - x (30) 

Finally we recall the expression of the inner product between two one-particle states in 
Fock space as expressed in terms of wave functions. Let the one-particle wave packets \ci) 
and |c 2 ), corresponding to momentum-space amplitudes ai(p, () and a 2 (p, C)> De represented 
by the wave functions ip\{x) and ip2{x) as in Eq. (J29l) . Then the equal-time inner product 
can be written in several equivalent forms: 



M*,t) = {o\Mx)\a) = J2j ^ (2 ^ ) 3 32E.( P ) e^P ' X "^ (P)ta(p, CK(p ' ' (29) 
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13|, 



(c 2 |ci) = rf3 P a 2(P>C)ai(P,C) = 5^ / rf3xS 2( x >C)ai(x,C) 

= y"d 3 x^ 2 (x,t) 7 > 1 (x,t) = y"rf 3 xV|(x,t)Vi(x,t). (31) 
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Notice that for a non-interacting field the value of this inner product is independent of time. 

III. THEORY OF DIRAC AND MAJORANA NEUTRINOS 

Given a four-component spinor ip(x) it is always possible to project out its left-handed 
(LH) and right-handed (RH) parts an d ^r{x) by means of the projection operators 

Pl = an d Pr = ^2~^ respectively Thus, 

Mx) = p L^) , M*) = p R^( x ) ■ ( 32 ) 

In the following we indicate a LH field for a neutrino of flavor a as Ul >oc (x). It is this field 
that enters the Lagrangian of the Standard Model. 

For Dirac neutrino fields, a mass term cannot be described by using the v^ a {x) field 
component only but a RH field v^ a {x) is needed as well. The term describing Dirac neutrino 
masses is [6] 

= ^2 v R,<* M oP u l,I3 + hx - ' ( 33 ) 

a(3 

where M a p is the mass matrix. 

In order to diagonalize M a/3 , one introduces two unitary matrices U and V jf], Q, in 
terms of which 

M a(3 = VamiU;, . (34) 

The quantity is the mass of the z-th eigenstate. The flavor eigenstates are rotated into 
the neutrino mass eigenstates and vrj as 

VL,i = U ai U L,a , VR,i = ^ K>«,a • (35) 
a a 

The mass term becomes 

£ m = ^ m i v R,i u L,i + h.c. = ^2 m i v i u i ■ (36) 

i i 

Here Vi = vr^ + ul,i- For future application we use the following parametrization for the 
two flavor generation mixing matrix U 

cos 8 sin 8 , 

U=\ I (37) 

— sm 8 cos 8 
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For Majorana neutrinos, it is possible to construct a mass term in the Lagrangian by 
using just one chirality, say LH, for the field v. The most general Majorana mass term is 

C = \ E ("LrfMapVLj, + h.C. . (38) 

Here {vl) c = CV[ T is the charge-conjugate of vl, and it is a RH field. In this case, in order 
to diagonalize the mass term in Eq. ( 1381) . only one unitary matrix W a i is needed, 

M aP = Y,W: a m l W; i . (39) 

i 

The rotated basis describing the left-handed mass eigenstates is 

"z,< = E ( 4 °) 

a 

In the new basis, the Lagrangian term in Eq. ( |38i) becomes 

C = ^E mi Ki) C ^,i + h.c. = ^m < x < Xi. (41) 

i i 

Here Xi — v l,% + (^z,*) is a Majorana field. 

The mixing matrix W of N Majorana neutrinos contains N — 1 physical phases in addition 
to the angles in the mixing matrix U of the Dirac theory (see e.g. [6(). This is due to the fact 
that the Majorana mass term in Eq. (1381) is not invariant under the U(l) global symmetry 



Following 



15l | it is possible to write the Majorana mixing matrix W as the product 
W = UD of a unitary matrix U similar in form to the Dirac mixing matrix in Eq. fl37|) and 
a diagonal unitary matrix D = diag(l, e 4 ^ 2 , e^ 3 , . . .): 

W ai = U m e^. (42) 

For the two-flavor case, one has explicitly 

(cos 6 sinOe 1 ^ \ 
■* ■ (43) 
— sin 8 cos 6e l<p I 

Finally, we recall what happens to the weak interaction vertices because of the change 
of basis from interaction to mass eigenstates. We take as an example the lagrangian for the 
weak interaction of neutrinos, charged leptons and W gauge bosons 

Cint = gi'L,a'flL,aW lt + h.C. (44) 
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FIG. 1: The Feynman vertex describing the weak interaction between a charged lepton l a , a charged 
vector boson Wn and a neutrino mass eigenstate of mass m,. The formula next to the diagram 
gives the strength of the coupling and depends on the mixing matrix U a i. 

In the last expression, Ul j0i is a LH neutrino field, which can be Dirac or Majorana, Il j01 is 
the LH spinor field for the charged lepton of flavor a, is the charged gauge boson field 
and g is the weak coupling strength. In the mass eigenstate basis, this lagrangian becomes 

Ant = gvL,iU a i'flL, a W lt + h.c. . (45) 
The vertex coupling a lepton ix,a of flavor a to a neutrino Vi of a given mass can be read off 

_____ _ 5 ^ 

Eq. (j45p as —igJJoa^—^-. This vertex is shown in Fig. [TJ 

IV. OSCILLATION FORMULA FOR DIRAC NEUTRINOS 

We now consider the following experiment, which resembles the actual operation of long- 
baseline neutrino experiments using accelerator beams. Imagine a burst of neutrinos created 
with a specific flavor a by some mechanism of production. We are not interested in the 
characteristics of the production process because, following standard quantum mechanics, 
once the initial conditions are given the propagation of the wave packet is determined. The 
source is fixed in the laboratory frame of coordinate Xp and the burst occurs at time tp 
in this frame. We set a detector at the point and we gather information about the 
neutrino flavor (3 at time tp> = tp + At. If desired, we could turn on the detector for a 
certain time AT^ e4 , in which case w would integrate the oscillation probability over the time 
interval AT^ ef . Again we do not take into account any theory concerning the mechanism 
of detection, but only assume we want the probability that the neutrino will be in a given 
final state at time tp>. Because of the phenomenon of neutrino oscillations, the flavor (3 of 
the neutrino burst at detection can differ from the flavor at production a. We are interested 
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in computing the probability to observe the flavor (3 at detection if the neutrino had a 
specific flavor a at production. Our final formula depends on the time delay At and on the 
characteristics of the mass eigenstate wave functions in momentum space Oj(p); it will be 
valid in general for any shape of the wave functions. 

In Section [VI] we specialize the wave functions Oj(p) assuming Gaussian distributions with 
average momenta P, and dispersions a p . The dispersion in momentum space is related to the 
uncertainty in the localization of the source (or detection point) by the relation a p a x = 1/2. 

We now proceed to derive the neutrino oscillation probability Vp a using a complete 
quantum field theory treatment. Following Eq ( 1261) we take the wave function describing 
the mass eigenstate of mass rrii at production as 

v - Ax) = £ / vAm^ CpMp - CpK " ■ (46) 

Here the extra label (P) indicates the production state. The wave function describing the 
flavor a is the linear combination of the wave functions Eq. (I46p 



^£4^,p(x). (47) 



We use the Dirac propagator Si(x — y) to evolve the wave function Eq. (j4"6"j) from the 
production time tp to the detection time to = tp + At. The wave function Eq. (146]) evolves 

as 

ViA x ) = i[d 3 ySi(x - yh\ P (y). (48) 
The propagator Si(x — y) is given by [11] 



(49) 



As in the laboratory frame the time delay is At, we rewrite Eq. (149]) as 
We integrate over ko using 

e -ik At e -iEi(k)At 

dkn— — = — 27u — - - . . (51) 

where we assumed that At > so the integration contour encloses the region Imfco > 0. 
This accounts for the propagation of particles to be forward in time. This gives 

Q.(AT X - V) - -7 ! ^ ^U ffl ' -»E,(t)Ai+ik.(x-y) /cnN 

5,(Ai,x y;_ i J (2?r)3 e . (52j 
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£ J ^ApXp)^ mSMr e-' g 'W" +A " + '-- (53) 



Notice that fcoAi i n Eq. (!50|) has now become Ei(k)At. The evolved wave function follows 

as 

u i>P (t + At, x) = % J d 3 ySi(At, x - y)j%, P {^ y) = 
= Y [ d 3 kd 3 pa lP (p,(p) ^ht±J^T ^M e -*g.(fc)At-»g,( P )wk.x f d 3 ye Hy(k-p) = 

(7% + m^)7°M(p, Cp) o -iB i (p)(t + At)+ip.x 

2^(p) v /(2vr) 3 2E i (p) ' 

In the last step we used J d 3 y e _4y ( k_p ) = (27r) 3 5(k — p). Now the following spinor relations 
hold 

7%7VP, Cp) = [-7°m + 2£7f(p)] u(p, Cp), (54) 

and 

(7>m + mh°u(p, Cp) = 2^(p)«(p, Cp)- (55) 
This gives the neutrino wave function after a time At as 

u hP (t + A*, x) = W /fn == , at,p(p, Cp)«(p, Cp)e-^W( t+A *) + ^- x , (56) 



Cp 



which is obtained from Eq. (j460 by letting t — > t + At. 

Notice from this detailed calculation that the propagator rightly contributes only the 
time dependent factor e -*^(p) A * ; which is the matrix element of the evolution operator in 
the energy representation. There is no extra term e *p ( x D- x p), The latter term will appear 
in the production and detection amplitudes when the respective wave packets are centered 
around the production and detection points. 

We now define the wave function centered around the detection point for mass eigenstate 

v . D ( x ) = y I d3p ai D ( p , CdHp, ( D ) e -^(p) i+J P- x . (57) 
The wave function describing the neutrino wave packet with flavor (3 at detection is 

V/3,d{ x ) = ) UpiVi,D{x)- (58) 
% 

The probability amplitude for the neutrino oscillations is defined as 
A af3 = J d 3 xu^ D (t,x)^°u at p(t + At,x) = ^2Up j U ai f d 3 xz/ j>D (t,x)7 z/ i>P (t + At,x). (59) 



12 



Inserting Eqs (|56|) and (1571) in Eq. (}59l) . the spinor normalization in Eq. ([6]) gives ( P = ( D = (, 
and we obtain 

-V = J2 U h U -i [ rf3 P <d(p, CKp(p, C)e"^ (p)At . (60) 

Eq. (160]) is our general result, obtained using basic principles of quantum mechanics only. 
No approximation was used to derive this equation, which is valid in general for any form 
of the wave packets ai tP (p, £). The production and detection mechanisms fix the form of 
the amplitudes aj,p(p, C) an d ch,d(p,C), respectively. Eq. ( 1601) gives the right probability 
amplitude with the correct normalization when At = and P = D. In this case the 
integrand reduces to \a it p(p, ()\ 2 d 3 p which is the probability of finding a particle in the 
infinitesimal volume of the momentum space d 3 p centered at p. 

We now discuss the properties of the neutrino oscillation probability Vp a . This is given 

by 



V, 



/3a 



\-Af3a\ 2 



J2 U k U ^ J d 3 pal D (pX)a hP (pX)e-^ ip)At 



(61) 



An interesting relation is obtained for the sum J^'P/fo, which represents the probability 
of observing a neutrino after time At independently of its flavor. Using Eq. ( 16T1) and the 
unitarity of the mixing matrix U, one finds 



(62) 



where V% is the probability of observing the z-th mass eigenstate if there were no mixing 
between neutrinos, 

2 



/ rf3 P a i,D(P>CKp(p,C)e" 



4Ei(p)At 



(63) 



Eq. ( |62|) states that the probability Y2p 'Ppa °f observing a neutrino in any flavor after time 
At equals the weighted average of the probabilities Vi of observing a neutrino in each mass 
eigenstate, weighted by the probability \U a i\ 2 of being in that eigenstate in the initial state. 

Eq. (1621) applies in general, for coherent and non-coherent cases, and for oscillations or 
lack of oscillations. As illustrated in Section IVIIII below, it allows a clear understanding of 
the phenomenon of decoherence, and of the question of the arrival time of the neutrino. 
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V. OSCILLATION FORMULA FOR MAJORANA NEUTRINOS 

Here we compute the probability amplitude for Majorana neutrinos and review the well- 
known result that there is no difference from the Dirac theory . The flavor-eigenstate 
wave packets are 

\c a , p ) = Y^ W ai\ c i,p) =$^W«i J d 3 P a p(pX P ) \p,( P ,i) , (64) 
i i,Cp ' 

where p = P,D labels the production and detection points and the states |p, C, i) are now 
defined in terms of the creation operators S|(p, Q of the Majorana theory. The flavor wave 
functions are 

XaA*) = YsW a MHx)\^ P ) = W ™J ^( 27r t )32^.( p ) ai ' p(Pi • (65) 

Following the same steps in the derivation of the oscillation formula Eq. fl60|) for Dirac 
neutrinos, we arrive to the following analogous expression in the Majorana case: 

Aap = J d 3 xxl i2 (x)xa,i( x ) = Yl W f* W <* J d 3 xxl,2(x)XiAx) = (66) 

= E W h W ™ [ d3 P<D(P, CKp(p, ()e~ lEdp)At . (67) 

When the matrix W is parametrized as in Eq. (142j) . one finds 



WZiWci = iU^YiU^) = UfrUa. (68) 

Thus 

A* = Y, U h U *i [ rf 3 pa* D (p,C)a l ,p(p,C)e-^ (p)A *. (69) 

This expression is identical to Eq. flBDl) for Dirac neutrinos. The physical phases fa play no 
role in the theory of neutrino oscillations and it is impossible to distinguish between Dirac 
and Majorana neutrinos on the basis of their oscillations phenomena (see e.g. 6|, [7J]). 

VI. GAUSSIAN WAVE PACKETS 

A minimal-dispersion Gaussian wave packet with average position X, average momentum 
P, spin polarization ( and spatial dispersion a x is represented at time t = by 

axW,0 = j^r^^^ . (™) 
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where 5^ is the Kronecker symbol which fixes the polarization along an axis £. Its 
momentum-space probability amplitude is a Gaussian with dispersion a p = l/{2a x ): 

«x,p,c>, = ^|pe-^e--- x . (71) 

Notice that the term e _jp X is a function of momentum, not position, and arises from dis- 
placing a wave packet whose average position is at the origin into one centered at X. 

We now specialize the production amplitude by assuming that the distribution in mo- 
mentum space aj i p(p, (p) is represented by a minimal-dispersion Gaussian wave packet with 
average momentum P^p, average position Xp, spin polarization along a fixed direction (p 
and momentum dispersion <r Pi p: 

a»,p(p, Cp) = ax,p,c,a p (p, 0- (72) 

A similar equation holds for Oj .p(p, Cd) by replacing P ^ D. The Fourier transforms of these 
packets are centered at the production (Xp) and the detection (Xp>) point, respectively. 



As discussed in 



BE 



9| Gaussian wave packets are just approximations of the true distri- 
butions, whose shapes will depend on the details of the production and detection processes. 

Using the distribution (1721) in momentum space we get the following probability amplitude 
for the process 

^-k^U h U m J ^Jj^ e . (73) 

The dependence on the production and detection points Xp and is now explicit because 
of the choice of the Gaussian distribution fl72|) . In Eq. fl73|) we introduced the notation 
AX = X D - Xp. 

Combining the quadratic terms in the exponential, one obtains the known result (see 
e.g. B]) that the relevant momentum dispersion o 2 vPD is a combination of the momentum 
dispersions at production and detection, 

^ = jr + jr- (74) 

°pPD °pD °pP 

The separation of the wave packets will be governed by the combined dispersion <j p pd- In 
particular, the spatial coherence of the source, as parametrized by a p p, enters the determi- 
nation of wave packet separation. 
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In the illustrative case in which the initial and final polarizations, the average momenta, 
and the momentum dispersions are the same, i.e. £d = (p, P^d = Pi,p = Pi and a v p, = 
a p p = Up, we have 

This is the basic expression we use below to illustrate neutrino oscillations and decoherence. 

We now expand the energy term Ei(p) appearing in Eq. fl75l) to second order in the 
momentum q = p — P,. We find 

«<p) = «<fl) + |^ + ^|«|' + 0(|S|'). (76) 

This expansion is a good approximation when -J-^py *C Ei(Pi). Since the momentum dis- 
tribution is Gaussian with dispersion a p , we have \q\ ^ u p and the expansion in Eq. (IT6]) 
applies for 

2 

777 

o p « P t + -±, (77) 

where Pi = |Pj|. The condition is valid for non-relativistic neutrinos (for which o v <C m; 
but <jp can be smaller or larger than Pi), while in the ultra-relativistic limit it is equivalent 
to a v <C Pi. This conclusion differs from that in [7], where this approximation is said to be 
valid only in the ultra-relativistic limit. 

The first order term in Eq. f!76|) is responsible for separation of the wave packets, while 
the second term is related to their spreading. To see this, we integrate Eq. (1751) with the 
expansion ( |76l) up to second order and find 

A/3a = ^r7 Q ^e-^^ A '- AX ) 2 ^^) AWP '- Ax . (78) 
i a v 

Here 

v - = m> (79) 

is the group velocity of the z-th mass eigenstate, and 

= YT^m [1 " Mt)] ' (80) 

where 
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The quantity c p (t) is the time-dependent dispersion of the wave packet in momentum space, 
in analogy with the spreading of a Schrodinger-like wave packet [161 ] . 

At small At <C [mf + P 2 ) 3 / 2 / (<7 2 m 2 ) (for all i), the spreading of the wave packets can be 
neglected, <J p (i) ~ a v . Then the oscillation amplitude and the detection probability become 

Apa = £ UfrUate-^^^^*^-**, (82) 

i 

V 0a = J2 ^^/^'^^ (83) 

Here cr x = l/2a p is the dispersion in position space. In these formulas the product of the 
IPs and the oscillating exponential are the usual expressions for plane waves, and the last 
exponential is the envelope describing the superposition of the two Gaussian wave packets. 
In Eqs. ( 1781 and ([52]) . each mass eigenstate i contributes a phase 

= £i(Pi)At - P, ■ AX. (84) 

We want to point out that the quantity At appearing here is not the time of flight of the 
neutrino mass eigenstate i. Rather it is the time delay between production and observation 
of the neutrino packet. This time delay is fixed by the experimental situation, and is 
independent of the neutrino mass eigenstate. On the contrary, the time of flight for the i-th 
mass eigenstate is given by Ati = AX/vi, and is different for different masses. If one were 
incorrectly to replace At with Ati, one would obtain a spurious factor of 2 in the oscillation 
formula (see e.g. |7j, Section 8.4.2, and references therein). Our expression does not contain 
such a spurious factor of two. Moreover, it is relativistically invariant (if the spreading of 
the wave packet is neglected). 

VII. TWO NEUTRINO FLAVORS 

In the following we illustrate the probability in Eq. fl83]) by considering two neutrino 
flavors only, e.g. \x and r flavors, neglecting the mixing with the electronic flavor or with 
other hypothetical sterile neutrinos. Imagine a muon neutrino produced by some mechanism 
at a source located at X = Xp. This muon neutrino is described as a superposition of two 
neutrino mass eigenstates wave packets of mass m x and m 2 with average momenta Pi and 
P2, respectively, starting at the same time t — with the same average position X. In this 
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scenario we can compute from Eq. (I83I) the probability V^ T of observing a tau neutrino 
after a time At with a detector placed at distance AX from the production point. Using 
the parametrization in Eq. ( |37l) . the oscillation probability formula in Eq. (1831) specializes 
to 



7? = sin 2 (20) 



z 1 -z 2 \ , / AEAt - AP ■ AX" 

+ ziz 2 sin 



^AEAt — AP ■ AX\ 



(85) 



2 

As explained before, the same expression is obtained for both Dirac and Majorana neutrinos. 
In Eq. ( 185]) we have defined 

Zi = e -*(AX^ (g6) 

Physically, z 2 is the overlap integral between the initial and the final Gaussian wave packets 
for the i-th mass eigenstate when the spreading of the wave packet is neglected, as in Eq. (I83j) . 
Notice that — 1 at Ax = At = and on the classical trajectory AX = VjAt, and Zj < 1 
elsewhere. The product z\z 2 sin 2 (2#) is the overlap integral between states of the same flavor 
at production and detection. Oscillations are present only when the latter overlap integral 
is different from zero. 

The probability to observe the same flavor /i can be similarly obtained as 

V ^ = (z\ cos 2 9 + z\ sin 2 9) 2 - Zl z 2 sin 2 (2fl) sin 2 ^ AEAt - AP ■ AX ^ (g?) 

From the last two equations we see that the total probability of detection of any flavor, if 
we started with the flavor /i, is 

= V^ T + = z\ cos 2 9 + z\ sin 2 9. (88) 

This probability is one when Ax = At = and when computed along the classical trajectory 
AX = VjAt. Away from the classical trajectory, the probability of observing a neutrino is 
smaller than one. This may seem to be an unusual result, since one would expect that 
the probability of observing one of the two flavors has to be 100%. However, what we 
are computing here is not the probability that if we start with a neutrino at the origin 
we observe it to be somewhere at a later time, which correctly is one. Instead, we are 
computing the probability that an initial Gaussian neutrino wave function centered at the 
origin at time t = has become a final Gaussian of the same dispersion centered at AX 
at a later time At. This probability is given by the overlap integrals zf between the initial 
and the final Gaussians. Eq. (IHHj) for the total probability of observing any flavor at a later 
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time is then easily interpreted as the sum of the probability of producing and detecting 

each mass eigenstate, eigenstate 1 being produced with probability cos 2 9 and observed 

with probability zf, eigenstate 2 being produced with probability sin 2 ^ and observed with 

probability z\. This is a special case of the general result expressed by Eq. ( 1621) . 

In the plane wave limit, o v — ► 0, o x — > oo, and Zj — > 1. Eqs. ( 18514881) reduce to 

r, / AEAt — AP • AX\ , . 

P^ r = sin 2 (20) sin 2 ( J , (89) 

9 / AEAt — AP • AX\ , s 

7V„ = 1 - sin 2 (20) sin 2 f J , (90) 

V, = 1. (91) 

Eqs. ( I89H9T!) contain the relativistically invariant combination AEAt— AP-AX = AP^AX^, 
where four-vector notation has been introduced. 

In the ultra-relativistic limit Pi ~ E — m 2 / (2P) where E pa E\ pa E-x to the first order. 
In this limit, AE = 0, AP = Am 2 /(2P), and 

P^ T = sin 2 (20) sin 2 f 4g j (92) 

= 1 - sin 2 (20) sin 2 f — j . (93) 

In an alternative derivation of the ultra-relativistic limit, we could set At ~ AX and use 
AP pa AP + Am 2 /(2P), which follows from P ps P; + m 2 i /(2E i ). We obtain AP^AX^/2 Pa 
AE(At - AX)/2 + Am 2 AX/(4P) pa Am 2 AX/(4P). Again this has the correct factor of 4 
in the denominator. 

However we remark that Eqs. fl89H9Tl) hold without assuming AP = 0. In particular, 
one does not need to assume that the energy or the momentum of the different wave packets 
are equal. Indeed, such assumptions are not Lorentz invariant, and would be valid in a 
particular Lorentz frame only. Depending on the particular process in which neutrinos are 
produced, there will exist a frame in which it is possible to equate either the bulk momenta 
or the energies of the eigenstates, but this makes little importance in our treatment as we 
already chose our reference frame by setting the relative distance of source and detector AX 
and the time delay At. 

It is convenient to define a coherence length between two wave packets. Their average 



position at time tp + At is given by Xj = Xp + VjAt, where = P*/ sj m? + P 2 is the wave 
packet group velocity. We define the coherence length L coh as the distance | (X 1 +X 2 )/2— X P | 
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from the source at which the two wave packets are separated by an amount equal to twice 
their spatial dispersion 2a x , that is |vi — v 2 | At = 2o x . The coherence length follows easily 

as 

T l v i+ v 2| {aA , 

Lcoh = i rcr x - (94) 

|vi - v 2 | 

The coherence length sets the distance from the source within which it is possible to observe 
oscillations, since it is gives the value of the overlap integral z 1 z 2 sin 2 (26*) (see also Figure (@J 
below). As pointed out in [l^] the coherence length is proportional to a x , the spatial disper- 
sion at production and at detection. It is thus related to the characteristics of productions 
and the resolution of the detector. As we stress in this paper, we are not interested in these 
details as once they are set the properties of the oscillation probability are determined. 



VIII. NUMERICAL ILLUSTRATION - ULTRA- RELATIVISTIC NEUTRINOS 

In this section we show the time behavior of the oscillation probability Vp a by means of 
a numerical illustration, limiting the treatment to the case of two flavor generations. We 
consider atmospheric neutrinos created by the interaction of cosmic rays with the Earth 
atmosphere. We take the following values for the parameters in Eqs. ( I85H871 ) . 

AX = 6000 km, Am 2 = 2.7 x 10~ 3 eV 2 /c 4 , E = 2 GeV, 6 = 45°, a x = 10 km. (95) 

We imposed the same energy E for both neutrinos. Since E 3> m«, the ultra-relativistic limit 
applies. In this limit one has Vi ~ 1 — m 2 / (2E 2 ), and the coherence length L co ^ becomes 

AE 2 

Lcoh ^ ~^2 ax ( for E > mi > 2 )- ( 96 ) 

For the values in Eqs. (195]) . one finds L co h ~ 6 x 10 22 km ~ 2 Gpc, comparable to the Hubble 
radius. So, oscillations are coherent on the scale of the Earth. 

For AX -C £ C oh, one can set z 2 = Z\ = 2, and the overlap integral z 2 factorizes in the 
oscillation formulas, 

^ 9 9/ ^ 9 f Am 2 AX\ , , 

V^ T = z 2 sin 2 (26) sin 2 — , (97) 



V = z 2 



1 - sin 2 (26) sin 2 



AE 
Am 2 AX 



(98) 



AE 

An important feature of the oscillation phenomenon in the ultra-relativistic limit is its 
dependence on Am 2 only. This fact makes an oscillation experiment insensitive to the 
individual values of the masses mi and m 2 . 



20 




5 10 15 

Distance traveled ( 1 OOOKm) 



FIG. 2: An example of the probability of detecting a neutrino in a Gaussian state when produced 
in a Gaussian state, plotted as a function of the distance from the source AX at a specific time 
At after production. We assume two neutrino generations with masses and mixing angles given 
in Eq. [95] and take At = 20 ms. The dashed line shows the probability of detecting either flavor: 
this probability peaks at AX = cAt = 6000 km corresponding to the classical trajectory for ultra- 
relativistic neutrinos. The two intertwined solid lines show the probability of detecting one or the 
other neutrino flavor; oscillations are clearly visible. 



From Eqs. (197H981) . the oscillation length is 



An 



E 



Am? 



(99) 



(This is the distance between two consecutive peaks in the probability of observing a cer- 
tain flavor.) In this relativistic limit, the ratio of the coherence and oscillation lengths is 
independent of Am 2 : 



J cdh 



E 



L osc 8ira p 



for E 3> m% ; 



(100) 



This ratio is just the relative momentum dispersion o v jE that determines the number of 
oscillation cycles before coherence is lost and the two particles wave packets separate in 
space. 

In Fig. El we plot V^ T and in Eqs. ( I9TH981) as functions of the distance AX at 

time At = 20 ms, for the parameters given in Eq. [95j The two intertwined solid lines show 
the probabilities of observing one or the other neutrino flavor. The maxima are spaced 
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by a length equal to L osc = 1.86 x 10 3 km, as per Eq. (}99l) . The dashed line shows the 
probability V^ T + of detecting either flavor. This total probability is the overlap 

integral z 2 and it peaks on the classical trajectory AX = cAt = 6000 km. As the neutrinos 
are ultra-relativistic, the probabilities of detecting different mass eigenstates overlap and it 
is not possible to distinguish between the two packets in this limit. 



IX. NUMERICAL ILLUSTRATION - DECOHERENCE 

We now discuss the phenomenon of neutrino decoherence in the case where the distance 
between the source and the detector is AX > L C oh- F° r this purpose we focus again on 
a two-generation model, with neutrino flavors /i and r, but we choose an admittedly non- 
realistic example of non-relativistic neutrinos created at the source with the same average 
momenta Pi = P2 = P. The parameters are as follows 

mi = 2 • 10~ 3 eV/c 2 , m 2 = 0.022eV/c 2 , P = 2.2 • 10" 4 eV/c, 

AX = 6000km, 9 = 45°, a x = 10km. (101) 



The oscillation probability is described by Eqs. (1851487]) . 

In the non-relativistic limit, the coherence length from Eq. (1941) becomes 

2(m 1 + m 2 ) 2 



J coh 



Am 2 

while the oscillation length is 



-a x , (102) 



mi + m 2 

2tt— — — . (103) 
Am z 



We have chosen the parameters in Eq. (110 II) so that L co h ~ AX. 

As can be seen in Fig. [3X oscillations are suppressed at late times, when the two wave 
packets do no longer overlap. The solid lines represent V^-> T and V^^, the transition 
probabilities of observing neutrinos of flavor r or fi, respectively, if the initial flavor is \i. 
These probabilities are bounded by 



8 -^l( Zl - Z 2 ) 2 < V^r < S ^(Zl + Z 2 f (104) 

**m(z x - Z2 f < < 8 ^(z 1 + z 2 y. (105) 



As v\ > v%, the probabilities of detecting a neutrino of flavor /iorr drops as \ sin 2 (26*)2; 2 
for At> AX/( Vl -v 2 ). 
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FIG. 3: The neutrino oscillation probability Eqs. (j85[) . (|87p as a function of the time delay At 
between neutrino production and detection, in the case of two-neutrino generations and for a 
coherence length L co h comparable to the distance AX between the production and detection points. 
The parameters used for the plot are listed in Eq. (jlOip . The dashed (blue) and dotted (red) lines 
represent the probabilities to observe the mass eigenstates of lighter or heavier masses mi or 
77i2, respectively. Notice that the lighter neutrino arrives before the heavier one. The solid black 
lines show the probability that the initial flavor /j is detected as r or /j after a time At. Oscillations 
occur only when the lighter and heavier neutrino wave packets overlap. 

The separation of the wave packets is clear in Fig. [3j If there were no mixing, since the 
two packets have the same bulk momentum P, the lighter wave packet (dashed blue bell-like 
curve on the left, peaking at ~ 0.2 s) would arrive sooner than the heavier one (dotted red 
bell-like curve on the right, peaking at ~ 2 s). With mixing, the heights of the peaks are 
suppressed in the manner of Eq. ( 1621) : the light-neutrino peak is suppressed by a factor 
cos 2 9 and the heavy-neutrino peak by a factor sin 2 9. Oscillations can occur only when the 
packets have a substantial overlap (out to ~ 0.8 s in our example). A detector with timing 
capabilities placed at AX > L coh would be able to distinguish the two mass eigenstates 
based on their different times of flight from the production point. The ability to discern 
the two mass eigenstates destroys the interference pattern, in a way similar to the famous 
double slit experiment. 

Another way of seeing the separation of the wave packets is to plot the detection proba- 
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bility as a function of the distance from the source at different times from production. We 
do this in Figure H] for three different times. To produce a clear illustration, we have chosen 
the following parameters: 

mi = O.leV/c 2 , m 2 = O.lleV/c 2 , P = l(r 3 eV/c, 6 = 45°, a x = 10km, (106) 

and At = 10 s, 25 s, and 50 s. The first snapshot at At = 10 s shows the two wave 
packets (solid blue and red lines) overlapping. As a consequence the probability of detecting 
a neutrino of any flavor is relatively high, as shown by the dashed line representing the 
envelope of the oscillation probabilities (not shown). The second snapshot at At = 25 s 
shows the wave packets when they are separated by a spatial distance equal to a x , which 
is half of the coherence length as defined by us. In the third snapshot, at At = 50 s, 
the wave packets are no longer overlapping; the probabilities of detecting a neutrino flavor 
mirror those of the mass eigenstates, but are suppressed by the probability of producing the 
specific flavor in the first place. 




Distance travelled (Km) 

FIG. 4: Illustration of how two wave packets produced at the same point of space-time with different 
group velocities gradually diminish their overlap losing coherence. Three snapshot at successive 
times are shown (from left to right): At = 10s, 25s and 50s. The solid lines show the probability 
of observing one or the other of the mass eigenstates. Notice how the peaks of the probabilities 
separate. The dashed lines show the envelopes of the flavor oscillation probabilities (not shown 
because they vary on a very short distance and would blur the figure if shown in their entirety). 
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X. DISCUSSION AND CONCLUSIONS 

In this paper we have derived the probability of detecting a neutrino with a specific 
flavor (5 by using basic rules of standard quantum mechanics. In doing this we only used 
properties belonging to the neutrino mass eigenstates i/j, the flavor eigenstates being just a 
linear combination of Ui according to the mixing matrix U (if Dirac) or W (if Majorana). We 
gave physical reasons about the choice of using the mass basis instead of the flavor basis, as 
the former diagonalizes the mass matrix in the hamiltonian. For the probability amplitude 
in position space we used the Fourier transform of the momentum probability a(p), as in 
the prescription by Newton and Wigner. This allowed us to find the correct expression for 
the oscillation amplitude, Eq.(l75i). 

We commented on the equal energy prescription often found in the literature on neutrino 
oscillations. Thanks to Eqs. (185|) and (157)1 . which express the oscillation probability in a 
Lorentz invariant way, we showed that this requirement is not necessary and oscillations are 
generally attained also when the energies and the momenta of the mass eigenstates are not 
equal. In the two- generation case, it may be possible to choose a Lorentz frame in which the 
neutrino energies are equal. In the general case of more than two flavors, this may not always 
be possible, depending on the details of neutrino production. If the production process is 



stationary, equal energies can be assumed (see 17|), but in general it is not assured that a 
specific frame can be found where all of the energies are equal. 

We have also stressed that the oscillation amplitude depends on the details of the pro- 
duction and detection processes only through the initial conditions for the flavor eigenstate 
wave packets. To illustrate the physics of the neutrino oscillations, we have shown some 
numerical examples for a two-dimensional flavor space. In particular, we have presented a 
simple explanation of the separation of the wave packets and the consequent loss of coherence 
of the oscillations. 
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